Effective field theory approach to fluctuation-induced forces between colloids at an interface 
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We discuss an effective field theory (EFT) approach to the computation of fluctuation-induced interactions 
between particles bound to a thermally fluctuating fluid surface controlled by surface tension. By describing 
particles as points, EFT avoids computing functional integrals subject to difficult constraints. Still, all infor- 
mation pertaining to particle size and shape is systematically restored by amending the surface Hamiltonian 
with a derivative expansion. The free energy is obtained as a cumulant expansion, for which straightforward 
techniques exist. We derive a complete description for rigid axisymmetric objects, which allows us to develop 
a full asymptotic expansion — in powers of the inverse distance — for the pair interaction. We also demonstrate 
by a few examples the efficiency with which multibody interactions can be computed. Moreover, although the 
main advantage of the EFT approach lies in explicit computation, we discuss how one can infer certain features 
of cases involving flexible or anisotropic objects. The EFT description also permits a systematic computation 
of ground-state surface-mediated interactions, which we illustrate with a few examples. 

PACS numbers: 05.40.-a, 03.50.-z, 68.05.-n 



I. INTRODUCTION 

Be it the electromagnetic field, the curvature of a surface, 
or the composition of an inhomogeneous mixture, boundaries 
or objects interacting with a field generally place constraints 
on the fluctuations of it. This is known to lead to interactions 
between the constraining objects. The earliest and most fa- 
mous example of this effect was discussed in 1948 by Casimir, 
who showed that two charge-neutral conducting plates in vac- 
uum attract due to the constraints they impose on the quantum 
fluctuations of the electromagnetic field |38 40 1. Today it has 
become customary to name such interactions after Casimir, 
even in cases where the underlying fluctuations are not quan- 
tum mechanical but thermal in origin I4ll443ll . In this arti- 
cle, we will investigate a particular incarnation of this effect, 
namely forces induced by thermal fluctuations between par- 
ticles bound to a fluid surface characterized by surface ten- 
sion HQ. 

The non-trivial aspects of such calculations tend to arise 
from the constraints that the extended objects impose on the 
partition sum. This issue is usually dealt with by pinning the 
field to the surface of the objects through delta-functions in the 
integration measure [47[|. A clear exposition of this method, 
applied to compact objects in fluid membranes and films, can 
be found in Ref. H48TI . This method was later refined in the 
context of the electromagnetic Casimir effect l5lll . In the lat- 
ter work, the constraints of the objects enter the interaction 
energy through their scattering matrix coefficients. In fact, 
usage of scattering approaches for electromagnetic Casimir 
interactions has a long history and has lead to some recent de- 
velopments l52l - l55ll in the field, concerning analytical results. 
Very recently, thermal Casimir interactions between biologi- 
cal membrane inclusions were studied in a similar spirit 15611 . 

In this paper we employ a different strategy to streamline 
the boundary condition issue, namely, reverting to a point par- 
ticle description. We emphasize right from the start that this 
is not an uncontrolled approximation for the shape of the ex- 
tended objects, because we can systematically rescue the com- 
plete finite size information by means of effective field theory 



(EFT). The key philosophy behind EFT is separation of scales, 
which is deeply rooted in the concept of renormalization (for a 
review, see Ref. 15711 ). Originally developed in the context of 
quantum field theory, it has been applied only quite recently 
in a purely classical framework by Goldberger et al. 15811 . who 
used it to study the gravity wave profile for inspiralling black 
holes. This particular incarnation of EFT has subsequently 
been utilized to derive not only new results in gravitational 
wave physics l59ll but also to calculate the leading order finite 
size correction to the Abraham-Dirac-Lorentz radiation reac- 
tion force law in classical electrodynamics 16011 . Both of these 
applications dealt with intricate boundary condition problems 
for classical non-fluctuating fields, whereas the present arti- 
cle extends the formalism to allow for finite tem pera tures. We 
recently introduced this generalization in Ref. ||5fJ; here we 
develop the formalism in greater detail, expand on our previ- 
ous results, and discuss further physical situations. 

Examples of point particle approaches employed to com- 
pute surface-mediated forces have appeared several times in 
the past l6lU63ll . in the context of biological membranes. 
However, in most of these references the final answer depends 
on the short wavelength cutoff of the continuum theory, which 
is ambiguous by construction. To alleviate the issue, the cut- 
off is usually fixed by reinterpreting it as a multiple of the 
inclusion size, but there are at least two problems with this: 

(i) The size of the objects embedded in the field and the short 
wavelength cutoff of the theory are generally unrelated, and 

(ii) even if the objects do have sizes of the same order as the 
cutoff, it is not clear how to treat particles with different radii. 
Instead, the proper way to free physical quantities from depen- 
dence on some ambiguous cutoff is via renormalization tech- 
niques. Our approach is constructed such that the cutoff never 
"contaminates" the interactions and therefore never inappro- 
priately restricts particle sizes. Another difference between 
earlier point-particle approaches and ours is the fact that our 
treatment is not approximate in that all finite size information 
is systematically retained in the point particle picture. Unless 
this is done, it would not even make sense to attempt calculat- 
ing higher order corrections to asymptotically correct results. 
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For the sake of clarity, we restrict in this paper to the case of 
rigid particles, e.g. colloids adsorbed at a fluid-fluid interface. 
In other words, within the region of the surface occupied by 
the particle, fluctuations are not only reduced, but are frozen, 
in a manner compatible with the boundary conditions at the 
circumference of the particle and its permissible rigid body 
motions. Furthermore, we will assume the shape of the parti- 
cles are such that they have a circular footprint on the surface. 
While these two assumptions allow us to carry out a more 
transparent discussion of the EFT formalism, we nonetheless 
outline what is entailed in relaxing these assumptions. We 
also briefly turn our attention to interactions that are not in- 
duced by fluctuations but by permanent shape-irregularities of 
the particles. The formalism allows a systematic approach to 
calculate these interactions without resorting to a superposi- 
tion assumption to disentangle the boundary conditions of the 
particles from one another. 

In the following section of our article, we begin our discus- 
sion by reviewing the nature of the surface and fluctuation- 
induced interactions we are interested in, as well as outlining 
our effective field theory approach. We then proceed in Sec- 
tion [III] with constructing the effective theory and computing 
interactions in an asymptotic expansion. First, we illustrate 
the computation by considering only one of the many terms in 
the expansion as an example in Section lTlIBI Afterwards, we 
derive the interaction free energy of a pair of particles as a full 
expansion, i.e. valid at arbitrary separation, and extract the 
leading near-contact behaviour. We then compute the lead- 
ing non-trivial triplet and quadruplet interactions, as well as 
briefly discuss the case of soft or non-axisymmetric particles. 
Finally, we show how the formalism can be used for problems 
where the particles impart a permanent deformation on the 
surface. To separate the main arguments from mundane tech- 
nicalities, we have collected details of certain calculations as 
well as some auxiliary discussions in four appendices. 



n. FORMALISM 

A. The surface 

We are interested in a fluid surface whose behaviour is de- 
termined by the surface tension energy a J dA, where the in- 
tegral encompasses the entire surface and a is the surface ten- 
sion. A nearly flat surface can be parameterized using the 
Monge gauge J64ll . describing the shape in terms of its or- 
thogonal deviation h from a flat xy base plane. This yields 
the energy functional 

n smi [h] = a j d 2 r^Jl + h?(r) , (1) 
s 

where we write the gradient with the partial derivative sub- 
script i. Note that the domain S of the integral is not the whole 
M 2 but excludes the regions defined by the projection of the 
colloids onto the base plane. The presence of the particles im- 
poses conditions on h at the boundaries of these regions. For 



concreteness, let us picture circular disks floating on the sur- 
face, which dictate the height of the fluid film anywhere along 
their circumference. 

Under the further assumption of a weakly-deformed sur- 
face, |V/i| <C 1, one finally arrives at the quadratic surface 
Hamiltonian 

n mzf [h] = ^a f d 2 r/i?(r), (2) 
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up to an irrelevant constant. The existence of this integral 
defines the function space h(r) belongs to. 

B. Fluctuation-induced interactions 

At a given temperature, the surface will fluctuate around its 
ground state shape, which minimizes the energy as defined in 
Eq. ©, while still subject to the boundary conditions imposed 
by the particles, i.e. the continuity of the surface along their 
circumference. These local constraints on the fluctuations will 
result in the partition function and the free energy being de- 
pendent on the spatial arrangement of the constrained regions, 
hence giving rise to forces. 

Formally, the free energy is given as 

/3F = - log J Vhe- fm ^ i[h] , (3) 

where the partition function is a functional integral over all 
permissible field configurations h{r), i.e. those compatible 
with the boundary conditions imposed by the particles. This 
integral is a Gaussian integral, on account of the Hamilto- 
nian (O being quadratic in the field, but its evaluation is not 
straightforward. Though one can formally write the free en- 
ergy © as ~ log dct K, this does not help with explicit com- 
putation, as the kernel IK = — V 2 (see Eqs. (f2]i andO])) is 
hard to diagonalize in this constrained function space Q. One 
wishes to free the functional integral in Eq. (O from the finite 
sized regions of constraint, such that the integral becomes a 
"straightforward" Gaussian. 

One possible means to this end involves extending the va- 
lidity of the Hamiltonian (f5]i to the entirety of R 2 at the ex- 
pense of modifying the integration measure by delta func- 
tions, appropriately chosen to make sure the functional inte- 
gral in Eq. (15]) only sifts out permissible field shapes h(r) liil - 
|4& |5ll |6l], [65D . In this article, we follow another path, which 
encodes the constraints in additional terms to the free sur- 
face Hamiltonian. This is different from the approach taken 
in Refs. ll65ll66ll in that it is not necessary to assume these ad- 
ditional terms to be small, which corresponds to the case when 



The integration kernel in Eq. {5) involves a Laplacian operator that acts 
on functions defined only over a subset S of IR 2 . Integrating by parts, 
Eq. (2) can be massaged into 2 d 2 rh(— V 2 )h + ^ Jq S Alhihdih 
which does not have the form i J d 2 rd 2 r' h(r)K(r,r')h(r') with a 
translationally invariant kernel K(r,r') = K(r — r') which could be 
diagonalized easily in momentum space. 
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the constraints on the fluctuation of the particles are not rigid. 
Even when one deals with rigid constraints, one can compute 
the free energy in a convergent perturbation expansion. This 
is achieved by means of an effective field theory (EFT). 



C. The effective theory 

As we have alluded to before, we want to describe the sys- 
tem with an effective Hamiltonian such that we can write the 
partition function as Z = J Vhe - ^""^ with the field vari- 
ations unconstrained. This allows the evaluation of the inte- 
gral over a Gaussian measure with the familiar translationally 
invariant harmonic kernel. To this end, we construct the ef- 
fective theory as H c s = H + AH where H is the free surface 
Hamiltonian 



H[h] = \a J & 2 



'h 2 dr) 



(4) 



and the perturbation AH is a collection of local terms specif- 
ically designed to capture the constraints, as we will discuss 
shortly. In this form, the effective theory describes a system of 
point-like particles embedded in an otherwise homogeneous 
surface. 

The free energy, measured with respect to the free energy 
of the unperturbed surface, is then calculated in a cumulant 
expansion 

@F = - log J Vhc-^e-P^ (5) 
= -log(e-^> = -X:^<H3AWr> c . (6) 
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The last step above can be taken as the definition of the cumu- 
lants denoted as (. . .) c , which involve many-point connected 
correlation funtions of the field. Since averages are computed 
in the Gaussian ensemble determined by the free Hamiltonian 
H, the many-point correlators break up into a product of two- 
point correlators 



1 



(h(r)h(r')} = —G(r,r') 
pa 



1 



47t/3ct 



log(r-r')', (7) 



due to Wick's theorem. Here, G(r, r') is the Green function 
or inverse of operator —V 2 , i.e. — V 2 G(r, r') = S(r — r'). 

The aim is that the effective theory reproduces the correct 
free energy; that of the original system with finite sized re- 
gions of constraint, or the full theory as we will henceforth 
refer to it. This is ensured by matching a set of physical ob- 
servables in the two theories, a procedure which will be ex- 
plored in detail further in Section llll Al 

In general, AH consists of all polynomials in the deriva- 
tives of the field h consistent with the symmetries of the prob- 
lem, evaluated at the position of each particle. As such, it 
has the form of an operator product expansion Il67ll involv- 
ing terms like A a h 2 (r a ), S /i? (r ), C a h 2 j (r a )h 2 (r a ,) etc., 
where the label a denotes a particle and r a its appropriately 



defined position, e.g. the center of a circular disk. The prefac- 
tors A, B, C, . . . are called Wilson coefficients and have to be 
fixed by the aforementioned matching procedure. An expan- 
sion of similar spirit was employed in Ref. l68ll for studying 
interactions involving spherical objects in a critical fluid. 

Let us briefly discuss the temperature dependence of the 
free energy expansion ©. Every cumulant, being the ther- 
mal average of a power of (3 AH, will involve many-point 
correlators that consist of a multiplicative combination of 
terms in AH. Recalling from Eq. (O that every pair of field 
occurrences in these many-point correlators carries a factor 
= k^T, each quadratic term in (3 AH can be counted as a 
factor each quartic term as and so on, in the term of 
expansion © they appear in. From this we see that, the tradi- 
tionally encountered form F = k^Tf (|r a |) of fluctuation- 
induced interactions ll44l - l48il6lll65ll66ll stems from quadratic 
terms in AH. This allows the interpretation of the interactions 
as those between induced "capillary charges," arising from 
thermal fluctuations around each particle and subsequent po- 
larizations due to the deformation that propagates from those. 
Similarly, one can see that higher-than-quadratic terms in AH 
give rise to an excess of factors of k-gT in the cumulants in 
Eq. ©■ Therefore, such terms in the effective theory pro- 
duce a free energy F ~ (k-B,T) p - 2 . These are higher order 
fluctuation corrections to the interaction. However, to treat 
these non-linear corrections properly, one has to also consis- 
tently relax the weak deformation assumption on the free sur- 
face Hamiltonian (|2j. This involves quartic (or higher order) 
terms in the free Hamiltonian that can also be treated perturba- 
tively. While such extensions do not affect the basic EFT for- 
malism, they require additional field theoretical sophistication 
that would distract from the basic idea we wish to communi- 
cate here, and hence we defer such non-linear corrections to 
future work. 



III. EFT OF FLAT DISKS ON A FILM 

After briefly outlining the main idea of our formalism, we 
will now proceed by applying it to a specific problem. We 
want to compute the fluctuation-induced interaction between 
colloids with a circular footprint on the base plane. For this, 
we write the effective theory as AH = J2 a AH a , where a 
labels the particles. For one particle located at r a , AH a is a 
derivative expansion of the form 



AH a 



(8) 



encoding a quadratic excess energy for deforming the field 
in the vicinity of that particle. Denoting the radius of disk 
a by R a , dimensional analysis shows that each Wilson co- 
efficient Ca^ scales as o-R 2e , such that the limit R a — > 
is well-defined. It is also worth pointing out that, applying a 
Hubbard-Stratonovich transformation to the partition function 
in Eq. ([5]), the effective Hamiltonian H c a = H + ^2 a AH a 
can be brought into the form of the effective action of Emig et 
al. JH 
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As we mentioned earlier, the terms in Eq. ([8]) represent in- 
duced charges localized at the points r a . Thus, the coefficients 
Ca will be referred to as polarizabilities. One can see this 
clearly by considering a generic term in Eq. dH) around an ap- 
plied or background field. Choosing r a = and dropping the 
label a for convenience, one substitutes h — > h hg + h into the 
term (1/2)C^ (<9 £ /i) and observes the term linear in h: 

C w d e h h z(0)d e h(0) = J d 2 rh(r) 

x cW/i bK (0) (-d) £ 6(r). (9) 



One can then identify the localized source 

J(r) = (-YcM&tfriO) 8*8(1 



(10) 



This source is an £ th order multipole induced by the £ th 
derivative of the background, created by thermal fluctuations 
or the field of another such charge. 

The constraints imposed on the motion of the particle (typi- 
cally by the experimental situation at hand) determine whether 
a certain multipole charge can be induced or not. For in- 
stance, if the particle is free to move vertically, there will be 
no monopole induced, or C' ' = 0, because the particle can 
accommodate a nonzero h(0) by moving up or down, without 
an excess energy cost. Similarly, the dipole polarizability C' 1 ' 
vanishes if there are no constraints on the tilt degree of free- 
dom of the particle, as it can then adjust freely to a nonzero 
dih(0). We will reprise these arguments in due time when we 
discuss the leading order interaction between the particles. 

The presence of a multipole like Eq. (TlOb on the surface 
creates a deformation around it, which is simply the following 
convolution of the charge with the Green function: 

5h(r) =-\ J d 2 r'J(r')G(r',r) 



-d e h hs (0)d e G(0,r) . 



(11) 



For example, Fig. Q] depicts the response of a horizontally 
fixed disk to a dipole background, h hg = r cos ip. The 
fluctuation-induced interaction can be viewed as the interac- 
tion between these induced multipoles mediated by the sur- 
face. 

Notice that the multipole responses of order £ > de- 
cay as r~ e , whereas the monopole response does not. In 
fact, the monopole response is not even square integrable, i.e. 
its energy Eq. (f2]i diverges logarithmically. However, this is 
not a property of the EFT point particle description, but of 
the surface Hamiltonian itself. Therefore, calculating inter- 
actions between objects with nontrivial monopole responses 
in K 2 , i.e. disks with pinned height, is intrinsically an ill- 
defined problem and requires regularization of the surface 
Hamiltonian — whether one is using EFT or not. The issue 
can be alleviated by considering a finite subset of R 2 or mod- 
ifying the theory with a damping term. However, since this 
would distract from the main idea we want to communicate in 
this article, we will defer this case to Appendix [D] and avoid 




Figure 1. The figure illustrates how a rigid, horizontally fixed par- 
ticle of radius R creates a local field deformation when exposed 
to a tilted planar background, h bs — r cos ip. The total defor- 
mation is the sum of the background and the induced deformation, 
8h — —(R 2 /r) cos ip. The discontinuity of the (normal) gradient of 
the total shape along the circumference of the particle is a dipole-like 
(~ cos ip) "capillary polarization" that can be thought of as sourcing 
the response Sh. 



monopoles in the main text by restricting to particles that can 
fluctuate up and down freely, i.e. particles for which C^ ' = 0. 

Having discussed the multipole responses, one more note 
for completeness is in order. Although one must build AH 
from all scalars that the symmetries of the problem allow, 



h 2 -- 



etc. are 



terms involving repeated indices, such as hi 
redundant. One way to see this is that the response Eq. ( fTTT i of 
such a charge to a background would necessarily be local, as 
Gu(0, r) = S(r), and hence such charges do not interact with 
anything away from them. 



A. Fixing the polarizabilities: matching 

Eq. ([8]) is the most general quadratic Hamiltonian compati- 
ble with the circular symmetry of one particle. The only thing 
we have not established yet are the exact values of the polar- 
izability coefficients, which we will attend to in this section. 
These values are fixed by matching a set of observables in 
the effective (point particle) and full (finite sized particle) the- 
ories. Although any suitably chosen set of observables will 
do, we will consider the most obvious one here: the polariza- 
tion response to an applied field. If one is dealing with more 
complicated physics and boundary conditions or shapes, a dif- 
ferent choice of observable may be more practical (e.g. for 
numerical matching) or even necessary. 

We gave the general expression for an £ th order induced 
multipole response in Eq. ( fTTT i. Clearly, the background that 
will induce this response is one that is £ th order in the coor- 
dinates, and the deformation decays like d l G(0, r) ~ r~ e . 
The polarizability should be fixed such that this EFT re- 
sponse is identical to the response in the full theory of the 
particle to the same background. Given the correspondence 
between our present Hamiltonian and classical electrostatics 
in 2D, the latter response is found analogously to the problem 
of a conductor in an external electric field. 
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If we write the derivatives in Eq. ( fTTT > explicitly, this ex- 
pression involves the index contraction /i ii s • (0)G,: 1 ...i e (0, r) 
between partial derivatives of the incident background and 
the Green function. The simplest choice of background is 
/i bg = ar e cos £ip. While a cartesian index contraction is not 
hard to do when the polarization is of order 1 or 2, a general 
expression is difficult to come by in this way for high orders. 
Luckily, if one carries out the contraction in complex coor- 
dinates (z, z) instead of (x, y), a general expression is found 
easily. The result is that the effective object will respond to 
a background of the form /i bg = ar e cos £ip by creating the 
deformation 



(2) (3) 



Sh(r) = -C (l) 



2 l £\(£-l)\acos£ip 



4tt<t 



(12) 



which is derived in Appendix lAl 

In the full theory, one has to solve the boundary value 
problem V 2 8h = for r > R with the boundary condi- 
tion h\ r=R = (h bs + Sh)\ r=R = 00. The full set of har- 
monic functions with square-integrable gradient has the form 
r~ l cos £ip with £ £ {1,2, . . .}. Hence the solution in the case 
of a background /i bg = ar e cos £ip is found to be 



Sh(r) = —aR 



21 cos £if 



(13) 



Comparison of Eqs. ( flZb and ( fT3l yields the complete set of 
polarizabilities 



4^i? 2 V 
2 l £\{£-\)\ 



1,2,3.. 



(14) 



for one object. For the sake of clarity, let us rewrite the com- 
plete AH by restoring the particle labels a: 

^w^EE^t^)] 2 , (is) 

a t=l 

where d l denotes d^d^ . . .di f . We now know the point- 
particle EFT Hamiltonian which — to quadratic order in the 
field — is rigorously equivalent to the full theory. 

As stated earlier, we do not intend to address higher-than- 
quadratic terms in the effective theory here. However, we 
would like to stress that the matching procedure is essentially 
the same even if such terms are considered. The steps taken 
since the beginning of Sec.|III]do not rely on the responses be- 
ing linear. What has been presented boils down to computing 
the one-point function (8h(r)) of the response 8h(r) of one 
isolated object to a background in the effective (point particle) 
and full theories. Generally, one matches such rt-point func- 
tions, or other observables that are functions or functionals of 



2 Strictly speaking, the boundary condition is h\ r _ R = a + b ■ r where a 
and b are free parameters describing the height and tilt (if applicable) fluc- 
tuations of the colloid. Their function is to satisfy the boundary condition 
readily without the need for an induced response for those backgrounds 
that the colloid can align with. 




(1) 

(1) 
(b) 



(1) 



(3) 

T3 -- 



(c) 



Figure 2. Diagrams depicting pair (a), triplet (b) and quadruplet (c) 
interactions stemming from the 4 th cumulant. Dashed lines repre- 
sent the world-lines of particles and have been labeled by the parti- 
cle. The number in parentheses at each vertex shows the multipole 
order at the associated point in the scattering process. Each propa- 
gator is subject to l\ + £2 partial derivatives where £\ and £2 are the 
multipole orders at the vertices it connects. Consequently, all these 
diagrams are interactions of order 2(1 + 2 + 3 + 1) = 14 in the 
inverse inter-particle separations. Also note that we show only one 
of three possible 4 body connections in (c), as our aim at this stage is 
just to work an example and not to be complete. 



them, across both theories to fix the Wilson coefficients. For 
the linear case, it turned out to be very straightforward to fix 
the coefficients independently of each other by matching the 
aforementioned one-point function. If nonlinearities are also 
accounted for, one may need to match more n-point functions 
to obtain a sufficient number of linearly independent matching 
conditions to fix the values of the coefficients and there may 
not be a simple choice of background that automatically de- 
couples them. However, despite the tedium introduced by the 
nonlinearities, the procedure is not intrinsically problematic. 



B. Fluctuation-induced interactions 

With the effective theory established, we can now turn our 
attention to evaluating the cumulant expansion (O explicitly. 
Taking powers of Eq. ( fl~5l l one observes that the q th cumu- 
lant in Eq. © involves 2g-point connected correlation func- 
tions. This can be represented as connected Feynman dia- 
grams with q vertices and 2 links coming out of each vertex 
(see Fig. [2] for instance), as one can see by Wick contrac- 
tion of the 2q point correlator. Every vertex contributes a 
factor of f3C^/2 while every link affords a two-point cor- 
relator (or propagator) as in Eq. (0, acted upon by deriva- 
tives stemming from Eq. (1151 . Recalling that C™> ~ a, one 
can see that there will be no factors of j3 or a left in any 
cumulant, i.e. the dimensionless free energy j3F is purely a 
function of the spatial configuration of particles. Each cor- 
relator (d e h(r a )d m h(ri,)) ~ d l d m G(r a , rt) brings a fac- 



tor \r„ 



-(l+m) 



Consequently, an interaction among 
q multipoles of orders £\ to £ q will scale with a total of 
—2(£i + £2 + . . . + £ q ) powers of the separations between 
them. Accordingly, one knows how many multipole orders 
and how many cumulants must be retained, in order to achieve 
any given level of accuracy of an asymptotic expansion of the 
free energy, expressed in inverse powers of inter-particle dis- 
tances. 
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It is instructive to consider one example term from the cu- 
mulant expansion, to make these ideas clearer. The fourth cu- 
mulant will serve well since it contains pair, three body, and 



four body interactions, while the expressions are still not very 
cluttered. Let us first write down the fourth cumulant directly 
from Eqs. © and (1151 : 



^ (4) : = "irE E b ^S Cdd 2 [^h(r b f [d^h(r c )] 2 [^h(r d )] 2 ) c . (16) 

a.b e a ,l b 

c,d e ot e d 



The summations are over the particles and the multipoles on 
each. In this section, we will discuss only a few terms embod- 
ied in these sums, chosen as instructive examples (shown in 
Fig.0. 

We will begin by considering the pair interaction between 
particles 1 and 2, depicted in Fig. [2ja). There are two factors 
one needs to compute: (i) the number of times the relevant 

product of polarizabilities, namely CiC^C^C^, occurs 
in the sum ( TTBT l. and (ii) the multiplicity of the Wick con- 
traction of the 8-point connected correlation function. The 
first of these factors is simply the multinomial coefficient 
(211) = 2!i!i! ' smce wna t we are looking for is the factor in 
front of a 2 bc in the expansion of (a + b + c + . . .) 4 . Next, the 
multiplicity of the Wick contraction, or the diagram Fig. EJa), 
is 2 4 due to the possibility of flipping the legs of each vertexQ 
Thus the pair interaction of Fig. a) can be written as 



(1)^^(2)^(3) 



ijk^-* jkl^lmnp^mnpi 5 V 1 '/ 

where the shorthand G ab = G(r a ,rb) has been introduced. 
For completeness, we also present the result after doing the 
index contraction (see Appendix [Aj> and substituting polariz- 
abilities from Eq. ( TBl i: 



/?F(Fig.|2ta)) - -6 



R 1 R 2 



12 



(18) 



where r ab = \r a - r b \. 

The triplet interaction depicted in Fig.[2fb) requires consid- 
eration of terms in Eq. ( fTBT ) which involve the polarizabilities 
C| and This combination occurs L ^ J 
times in Eq. ( TTSI ), which means, also recalling that the mul- 
tiplicity of the connection is 2 4 , that we obtain for the triplet 
interaction, 



^(Fig.Hb)) 



1/ 4 
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(19) 



3 It is generally 2 q with q being the cumulant order, except for q = 2, where 
it is 2. 



After index contraction and substituting the polarizabilities, 
we have 



/3F(Fig.Hb)) = +12 cos ^3 



d2 nlO r?2 
2 H l H 2 ^3 



7 7 
r 12 r 23 



(20) 



where is the angle between r\ — and r 2 — r 3 , or in 
this case, the exterior angle of the triangle formed by the parti- 
cles at vi- The quadruplet interaction of Fig.HJc) is similarly 
found as 



/3F(Fig.Hc)) 



~ 4 
4! VI, 1,1,1 
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n2 p4 n6 n>2 

- 12 005(3^+4^ ^777 
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(21) 



with an analogous definition of the angles. 

In this section we have tried to illustrate the Feynman rules 
of our expansion by specific examples. It should be clear 
that computing interactions is a straightforward exercise in 
counting powers and diagrams relevant for the desired order 
in the inverse separation. The Green function products are 
easy to evaluate using complex coordinates as discussed in 
Appendix lAl and they turn out to be proportional to a product 
of powers of the distances and the cosine of a combination of 
angles associated with the geometry of the particle configura- 
tion, determined by the multipole orders involved. 

The reader may have realized that there seem to be diver- 
gent pair and triplet interactions in Eq. ( fTSI l due to one or more 
self-links in the form d n G(0) (c.f. Fig.0. These divergences 
are all power-like (~ r~ n as r — > where n > 2) and carry no 
physical information, as there is no non-trivial renormaliza- 
tion group (RG) flow. The rigorous way to deal with them is 
to add counter-terms to the Hamiltonian, aiming to cancel the 
unphysical divergences from the free energy expansion. How- 
ever, since one can find all these counter-terms by renormal- 
ization techniques, and the sole purpose of them is to remove 
the divergences from the free energy expansion, the upshot of 
this systematic RG treatment is that the divergences can sim- 
ply be dismissed. We have included Appendix [C] where the 
removal of one such divergence by a counter-term is explic- 
itly illustrated. 



7 




Figure 3. Example diagrams with unphysical divergent self- 
energies. 



I. Pair interactions 

Although every multibody order is included in the cumu- 
lant expansion ((6), a systematic enumeration of all relevant 
diagrams lies beyond the scope of this article. However, if 
only two particles are involved, the solution can be worked 
out without too much trouble and one can write down the com- 
plete asymptotic expansion of the free energy. 

Let us first establish that only even-numbered cumulants 
are relevant for pair interactions. The reason is that there 
is no physical {i.e. finite) pair term in odd-numbered cumu- 
lants. The easiest way to see this is to consider a pair dia- 
gram, i.e. one with connections between two world-lines such 
as Fig. EJa), with an odd number of vertices distributed along 
them: There is no possible connection (Wick contraction) 
without at least one self-link. As we discussed right before 
this subsection, such contributions to the free energy are un- 
physical artifacts which can be removed rigorously. Hence, 
pair interactions exclusively stem from even-numbered cumu- 
lants. 

To write the full expansion of the pair free energy, one must 
be able to do two things at an arbitrary cumulant order 2s: (i) 
enumerate all contributing pair diagrams with 2s vertices and 
(ii) evaluate the diagrams, which involves doing the product 
of 2s propagators [such as in Eqs. (fTTI i, ( f20b or (ETT il with any 
collection of multipole orders at the vertices. Luckily, both 
of these can be done. The propagator product is calculated in 
AppendixlAl As for the enumeration of relevant diagrams, it is 
most suitable to again illustrate it with an example: The sixth 
cumulant contains the diagrams depicted in Fig.H] These are 



m 
w 



m 






Figure 4. Different "stitches" of the sixth cumulant pair interaction. 



all the connected pair diagrams with six vertices, free of un- 
physical self-energies. As one can see by starting from any 
one of the vertices and "stitching" the vertices in all possible 
ways, there are 3!2!/2 = 6 of these diagrams, which general- 
izes to g s := sl(s — l)!/2 for the (2s) cumulant (except that 
<7i = 1). With each given collection of multipole orders at 
the vertices, the contribution of each diagram to the cumulant, 
or the propagator product, will be different (not in terms of 
power of r but its prefactor) since the partial derivatives will 
generally be distributed differently on the propagators. While 
this might seem to complicate things, remember that there is a 
sum over all possible collections of multipole orders out front 
[c.f. Eq. (H6H . Thus, the effect of these g s different topolo- 
gies is to repeat each scattering process g s times in the sum. 
Hence one can consider only one of the g s possible diagrams 
and multiply the result by g s . 

Taking everything together, the entire pair interaction can 
be written as 



0F = - 



E 



i (2 



2, U 




(22) 



where J2{e } denotes a sum over all possible collections of 
multipoles at the vertices. Using the exact expressions for the 
polarizabilities, Eq. (fT~4T >. and propagator products, Eq. (IA1 lb . 
to evaluate the diagrams, the pair interaction can be expressed 
in the following asymptotic expansion: 



(3F 



E7 



E 



n 



+i 



iy.Rl 



2L Q D 2L C 
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r 2{L a 



(23) 

where i 2s +i = h , L = £ odd . I, and L c = £ cvon 1 1, . 
Assuming identical particle radii to reduce clutter, the first few 
terms of this expansion are 
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697 
IT 



12 



(24) 



Note that obtaining these numbers does not involve anything 
more complicated than elementary algebra. The interaction up 
to order r~ 70 has been plotted in Fig. [5] A cursory observation 
of the plot reveals that for (center- to-center) distances larger 
than r sa 4, the lowest order term in Eq. (l24l i suffices. How- 
ever, for closer separations, the number of higher order terms 
one should include increases rapidly. We demonstrate this in 
the inset to Fig. [5] where r 5 %(P) is defined as the separation 
where the series d24T i terminated at (D(r~ p ) achieves an accu- 
racy of 5 percent. As the dimensionless distance d = r/R — 2 
between the edges of the particles is decreased, the power 
where the series can be truncated and still achieve less than 
5 percent error appears to asymptotically obey P « 8.2 
From this we see, for instance, that if the surface-to-surface 



(2) 



(3) 



(1) (1) 
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Figure 5. Comparison to numerical results of Ref. [45] on the attrac- 
tive force f(r) = —dF/dr. The dotted and solid curves correspond 
respectively to Ref. 14511 and our results (truncated at order r~ 70 in 
the free energy). The dashed line represents the lowest order asymp- 
totic solution, i.e. the dipole-dipole interaction. The inset depicts 
how the power P of the highest order term included in the expansion 
and the distance at which a 5 percent accuracy is achieved depend on 
each other. 



separation is a tenth of the particle radius, i.e. d = 0.1, then 
one must retain all terms up to 0(r~ 82 ) in order to achieve 5 
percent accuracy. 

The first two terms of Eq. ( f24l > are the dipole-dipole 
and dipole-quadrupole interactions also computed earlier by 
Lehle and Oettel 14511 . The same publication also contains a 
quadrupole-quadrupole interaction of order r~ 8 , for the case 
of tiltable colloids, of magnitude —9. Since the expansion 
d24l i was developed under the assumption of frozen tilt degree 
of freedom, the magnitude of the r -8 term includes additional 
interactions and is therefore different. Namely, these addi- 
tional interactions are the quadrupole-octupole and dipole- 
dipole-dipole-dipole interactions (c.f. Fig. [6]) with respective 
magnitudes —6 and —1/2, such that the resulting total magni- 
tude of ther -8 termis —9—6—1/2 = —31/2. Had we instead 
assumed the colloids to be free to tilt, their dipole polarizabil- 
ity C^ 1 ' would vanish — hence so would diagrams Fig. |6|b) 
and[6]c) — and the expansion d24b would have started with the 
quadrupole-quadrupole term — 9(i?/r) 8 . 

Lehle and Oettel also compute a numerical solution to the 
interaction for intermediate and close distances, where the 
first few terms of the asymptotic expansion do not suffice l45ll . 
As a confirmation of our result, we show in Fig.[5]the interac- 
tion found in Ref. l45ll by performing the functional integrals 
numerically, with our expansion d23l l up to order r~ 70 . The 
agreement is excellent, except upon approaching the contact 
distance r = 2R, where our truncated series remains finite 
while the numerical solution appears to grow without limit. 
The exact solution indeed diverges as 



(a) (b) (c) 

Figure 6. (a) The quadrupole-quadrupole interaction — 9(R/r)~ 8 . 
(b) The dipole-octupole interaction — 6(R/r)~ a . (c) The dipole- 
dipole-dipole-dipole interaction — (l/2)(i?/r) -8 . 

where d = (r — 2R) /R is a scaled surface-to-surface distance. 
The dominant term has previously been predicted by combin- 
ing the exact result for two parallel lines with the Derjaguin 
approximation ll45l lr59Tk the subleading contributions cannot 
be extracted in this way but follow by expanding the exact 
analytical solution which one of us recently found by using 
conformal field theory techniques lf7(ill . 

Since our formalism by construction produces an asymp- 
totic expansion in the inverse separation, it is clearly not the 
preferred method when it comes to near-contact separations 
(see also inset in Fig. [5]). However, it turns out that even 
our large distance expansion (l23l l is in a sufficiently manage- 
able form to permit extraction of the leading order of d25l — 
including the prefactor — by some further analytical manipu- 
lation. We illustrate this in Appendix |B1 This of course con- 
stitutes an application of our expansion outside its region of 
highest analytical convenience, but it vividly illustrates that 
(i) our expansion can indeed be driven up to arbitrarily high 
order and that (ii) the resulting predictions are correct up to 
the point where the exact solution diverges. 

2. Beyond the pair interaction 

When there are more than two particles involved, the free 
energy of interaction generally does not decompose into a sum 
of pair interactions but contains multibody interactions like 
the ones discussed in Sec. lIIIBl As illustrated earlier, simple 
power counting reveals to what order in the separations a cer- 
tain multi-body interaction might contribute, and evaluating 
the interaction poses little difficulty. 

Let us now discuss the leading order triplet interaction. 
A dipole-dipole-dipole term from the third cumulant is the 
lowest order triplet interaction conceivable. It is easy to 
see by power counting that this interaction would scale as 
r v2 r 2z r z\ ■ However, the strength of this interaction turns 
out to be identically zero. This is due to the harmonicity of 
the Green function of the problem, and applies to any dia- 
gram in an odd-numbered cumulant (please see AppendixlAl. 
Therefore, with this choice of surface Hamiltonian, many in- 
teractions that would seem to exist based on scaling arguments 
actually disappear^. 



4 This is not the case for a membrane, where the surface energy stems from 
bending elasticity (Till . 
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(i) 



(1) 



Figure 7. (a) The fourth cumulant diagrams contributing to the lead- 
ing triplet interaction, (b) The fourth cumulant diagrams contributing 
to the leading quadruplet interaction. 



Since no triplet interaction exists in the third cumulant, one 



must look for the leading triplet term in the next. The lowest 
order interaction is between induced dipoles and stems from 
the diagrams depicted in Fig. [71 a). We have already computed 
a subleading correction to this as an example in an earlier sec- 
tion, that involved higher order multipoles. The leading order 
is similarly computed and found to be 



optrip _ 
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(26) 



which becomes, after evaluating the propagator product — 
either by direct computation or using Eq. (IA1 II ) — and sub- 
stituting the polarizabilities from Eq. ( fl4b . 
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(27) 



The fourth cumulant is also where the leading order quadru- 
plet interaction resides. This results from the diagrams shown 
in Fig.UJb) among four induced dipoles. It follows as, 
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(28) 



Depending on the desired level of accuracy, corrections to 
these interactions as well as further multiplets may be re- 
quired. In this case, one simply identifies the relevant dia- 
grams or interactions by power-counting and the computation 
that follows is fairly straightforward as shown. 



C. On anisotropic and flexible objects 

We have carried out our discussion under the assumption 
of rigid and circular particles. An obvious extension is to- 
wards particles with internal flexibility or non-axisymmetric 
shapes. Although a thorough treatment of these extensions 
can be quite involved, there are nevertheless certain simple 
but rigorous remarks we can make at this stage. 

Not much changes when the particles considered are flex- 
ible (or the constraints on their motion are not completely 
rigid). Since the boundary conditions at the circumference 
of the particles change accordingly, fixing the polarizabilities 
will require the solution of a different boundary value prob- 
lem on the full theory side of the matching procedure, but the 
form of the effective Hamiltonian is exactly the same. The 
subsequent change in the values of the polarizabilities will re- 
sult in different interaction strengths. However, this does not 
affect one's ability to identify features of the interaction, such 
as what the leading power of inverse separations in a given 



multi-body interaction is. 

If the elastic properties of the particles are ever so slightly 
different from the surface, then the free energy expansion ac- 
quires another smallness parameter intrinsic to the polarizabil- 
ities. This renders high cumulants in j3F negligible compared 
to the second; a so-called weak coupling [54;] approximation. 
In such a case, being able to obtain a closed form result is 
more likely since only the second cumulant has to be taken 
into account. This is one of the limits considered in Ref. lf56Tl . 
These authors restrict to two scatterings, analogously to trun- 
cating the cumulant expansion after the second, which is jus- 
tified by the weakness of the perturbations. 

In the case of an anisotropic particle, the response of the 
particle to incident backgrounds depends on its orientation, 
which is clearly not the case for the response given in Eq. ( fTTT ). 
One now needs to consider polarizability tensors to describe 
the anisotropic response of a non-circular particle. The first 
few terms in the derivative expansion of a particle read 



l ij C ijkl h kl + 



(29) 



Fixing the polarizability tensors will require solving more 
complicated boundary value problems, but the lowest order 
interaction will still stem from a pair of induced dipoles (ten- 
sor c£. ) and therefore scale as r~ 4 , no matter what the 
shape of the particle might be. Or if the particle is free to 
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tilt, then the polarizability of any term involving an hi will 
vanish, thereby making the leading asymptotic interaction a 
quadrupole-quadrupole at order r~ 8 . These are indeed in line 
with the findings of Ref. Il46ll . where explicit calculations of 
the leading order interactions between rigid ellipsoidal parti- 
cles were presented. 



IV. ON CORN FLAKES 

So far, we have focused on interactions initiated by thermal 
fluctuations and have thus restricted to particles that do not 
deform the film in its ground state. However, it is straightfor- 
ward to include such effects in our formalism. This kind of 
situation naturally arises from an irregular three phase contact 
line between a colloid and the fluid-fluid interface it is trapped 
at l72ll . Macroscopically, we see corn flakes perturbing the 
surface of a bowl of milk in a similar way. Again, one would 
expect (and can indeed observe) that once these deformations 
overlap, they give rise to aggregation of the flakes. 

These perturbations are treated as "permanent charges" in 
the effective theory, beside the induced charges arising from 
constraints. The effective theory incurs local linear terms in 
addition to the quadratic ones upon this modification. The 
easiest way to determine these linear terms is again through 
an argument involving the response of each particle to a back- 
ground. Consider a particle whose contact line height pro- 
file has a nontrivial multipole h cl (ip) = rj£ cos (£ip — tai) 
around its center with a phase angle a? . When a background 
h hg (r) = r)tR~ r cos (£(p — lat) is incident on this particle, 
the boundary condition (h m = h out )\ r _ R is automatically 
satisfied by the fact that h hg (R, ip) = h cl (ip) and hence no net 
response will be triggered. Another way to phrase this is that 
the field of the permanent charge and the induced response to 
this specific background eliminate each other. Therefore, if 
this "preferred background" is denoted by p^ l '(r), the local 
linear terms representing the permanent charges can be en- 
coded into the quadratic induced charge terms as a shift: 



(30) 



a.t 



Note that the values of the polarizability coefficients are the 
same, whether the preferred shape of the particles are flat or 
curved, as long as the constraints are the same. Notice that 
in the absence of (vertical) external forces, the vertical move- 
ment of a particle is free and hence = 0. Similarly, in the 
absence of (horizontal) external torques, the tilt motion is free 
and CW = 0. In the following, we will assume that this is the 
case and therefore the multipole expansion starts at £ = 2. 

From Eq. ( f30T > one sees that the linear terms in the effective 
Hamiltonian are 

-Y^CWd«pW(r a )d e h(r a ), (31) 

a,i 

which encode the permanent sources 

-Y,(-yC { a n d e \s(r-r a )dV\r a )} . (32) 



In the diagrammatic expansion, the linear terms correspond to 
vertices with only one link attached. Therefore, interactions 
involving permanent charges are open diagrams initiating and 
terminating at these. Such diagrams have one link less but the 
same number of vertices compared to the fluctuation-induced 
interactions, and power counting shows the free energy due 
to them will scale as F ~ aR 2 , i.e. with no temperature de- 
pendence. These ground state contributions together with the 
fluctuation part, ~ k-gT, of the previous section result in a free 
energy of the form F = ciq + a\T, which allows us to view 
this separation into ground state energy and fluctuation correc- 
tion as equivalent to a separation into the energy and entropy 
terms in F = E — TS. Observe that this equivalence would 
not hold had we considered anharmonicities in the theory and 
thus furnished the free energy with contributions higher than 
linear order in the temperature. 

For the sake of clarity, let us consider a pair of particles 
with a saddle shaped contact line deformation. That is, each 
particle will have a preferred background of the form p(r) ~ 
{\/2)kr 2 cos(2(/3 — 2a) around its center, where the angle a 
describes the orientation of the saddle. For a colloid under no 
external vertical force {e.g. gravity) or horizontal torque, this 
quadrupole deformation is the lowest order multipole contact 
line irregularity that can occur. The linear part of the effective 
theory follows as - C °- Pij(r a )hij (r a ), where 



cos 2a a 
sin 2a a 



sin 2a a 
cos 2a a 



(33) 



Notice the simple geometric interpretation of k a as the mag- 
nitudes of each principal curvature of the saddles. 

The leading interaction stems from the second cumulant as 
a direct interaction between the two permanent quadrupoles, 
Fig. [HI a), easily found as 



2 \ cf C< 2) 



2! VI, 1 



37rercos(2ai + 2a 2 ) 



P l A r i) G ljki Pki{r 2 ) 
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(34) 



where the pair was assumed to lie on the x axis to declut- 
ter the expression. The functional dependence of the inter- 
action implies both an attractive force and vertical torques 
on the saddle-like particles. The same potential was found 
by Stamou et al. 117211 . We have seen earlier that, between 
colloids that are free to tilt and fluctuate vertically, the lead- 
ing fluctuation-induced interaction is of order r~ 8 , which can 
only compete with the interaction (f34-b at close separations. 

Although we have restricted to a toy problem where the 
only permanent deformation is quadrupolar, it should be ob- 
vious how one would handle higher order multipoles of the 
contact line irregularity for accuracy at closer separations. 
However, one should note that the true interaction potential 
is not merely a superposition of direct interactions between 
permanent charges, similar to Fig. [8ja), but includes scatter- 
ings from induced charges, such as Fig. [8jb). The former is 
equivalent to assuming the total shape of the surface is given 
by a superposition of the deformations caused by each parti- 
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Figure 8. (a) The direct interaction E^ ~ r~ 4 between two per- 
manent quadrupoles. (b) The pair interaction E^ ~ r" between 
a permanent quadrupole and the quadrupole it induces on the other 
particle. This is the lowest order energetic pair interaction that in- 
volves scattering from an induced charge, (c) At the same order as 
diagram (b), there also exists a triplet interaction E^ involving two 
separate permanent charges and one induced charge. 



cle independently of every other, i.e. neglecting their bound- 
ary conditions, which is why it misses induced charges. Such 
a superposition approximation is captured by the second cu- 
mulant and holds up to the order the next cumulant begins to 
contribute. This is the order r~ 8 , where the third cumulant 
contributes the diagram in Fig. [Ub), evaluating to the repul- 
sion 

= lna{klR\ + kll%)?^ , (35) 

which is the lowest order pair interaction that involves a scat- 
tering from an induced charge. At the same order, one ob- 
serves the first multibody interaction as well, which is de- 
picted in Fig. |8lc), and evaluates to E^ = e[^ 3 + i%j + 

£312 where 

E [% = 9 wa ^&%h cos^a _ 2ai + 2a 2 ) (36) 

r 12 r 32 

and so on. If the particles possess higher multipole order con- 
tact line irregularities, merely taking these into account within 
a superposition approximation does not capture every possible 
correction systematically, since it misses contributions such as 
Figs. Ob) and[HJc) by construction. 

To summarize, although we primarily use the EFT formal- 
ism to facilitate the evaluation of a partition function, the 
ground state interactions mediated by the surface can easily 
be calculated as well. One can readily encode the permanent 
deformations imposed by the colloids as linear terms in the 
effective theory, upon which the ground state energy strictly 
separates from the fluctuation correction. Computations of 
this energy that rely upon the superposition of perturbations 
caused by the particles only hold for direct interactions of the 
second cumulant. 



V. CONCLUSION 

We have presented and discussed a formalism to compute 
the surface-mediated forces between rigid disks trapped at an 



interface, based on effective field theory. This formalism re- 
duces each finite sized particle to a point, which greatly sim- 
plifies computations. Still, these points are by explicit con- 
struction equipped with all the information necessary to cap- 
ture their finite sized counterpart's behavior, and hence the 
treatment is not approximate. In particular, their size is not 
simply recovered as the continuum theory's ad hoc ultraviolet 
cutoff — an approach which can neither conceptually satisfy 
nor be plausibly extended to more than one particle radii. 

The extension from the free surface theory to the effective 
theory is achieved by the addition of a derivative expansion 
AH, with extra terms localized at the centers of the particles, 
encoding the constraints they enforce on the fluctuations as 
induced capillary multipoles. Symmetries determine the form 
of this extension to the Hamiltonian, while one still has to 
fix the values of undetermined polarizability coefficients via 
a matching procedure. However, one gains insight into the 
interaction free energy even without fixing these coefficients, 
such as the leading order distance dependence of the force or 
the vanishing of certain many -body interactions due to the ker- 
nel of the theory. In this article we focused on circularly sym- 
metric particles, which allowed us to determine the effective 
Hamiltonian completely, i.e. to all derivative orders. However, 
in cases where this is not possible and a certain level of accu- 
racy — in powers of r — is desired, knowing only a fraction of 
all the polarizability coefficients will suffice. 

After the effective theory has been established, the interac- 
tion free energy is then computed in a cumulant expansion. 
iV-body interactions are identified as Feynman diagrams link- 
ing N world-lines and are computed in a fairly straightforward 
manner. The two expansions, the cumulant expansion and that 
in the local derivatives of the field at the inclusions, afford a 
transparent interpretation of every term that contributes to the 
interaction as a scattering process between various multipole 
moments of the induced charges, similar in spirit to Ref. 15111 . 
Unless both of these expansions are performed, consistent re- 
sults cannot be achieved. Also, owing to simplifications due 
to the harmonicity of the theory, the number of dimensions, 
and the symmetry of the particles, a general expression for 
any possible diagram was obtained in Appendix lAl We have 
used this expression to develop a full asymptotic expansion of 
the pair interaction. It was possible to partially sum up the 
expansion. We believe a complete summation of the series 
to obtain a closed form expression may be possible as well, 
though we have not succeeded here. 

As we have discussed, an extension to non-rigid and non- 
circular objects entails no major change in the formalism. 
The former only affects the strength of the polarizability co- 
efficients and therefore the interactions, while the latter re- 
quires consideration of tensor polarizabilities, and both fixing 
their values and computing interactions become a bit more in- 
volved. 

Moreover, we have presented results on ground state in- 
teractions between the particles as well as the fluctuation- 
induced interactions. This required us to add permanent 
charges to the effective theory, which is simply achieved by 
shifting the zero of the induced charge terms. The ground state 
interactions then strictly separate from the fluctuation correc- 
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tions and can be computed in the same manner we discussed 
for the fluctuation-induced forces. In addition to deriving the 
lowest order pair interaction, in agreement with known results, 
we illustrated how corrections — higher order or multibody — 
are computed. 

Note that the same formalism can be used to treat surfaces 
of different characters, i.e. different free surface Hamiltonians. 
A relevant example is the problem of surface-mediated inter- 
actions between inclusions on a biological membrane, which 
we will treat in a future publication 17111 . Additionally, fluc- 
tuation corrections of higher-than-linear order in U-qT can be 
systematically studied by increasing the anharmonicity of the 
theory, e.g. relaxing the weak deformation assumption on the 
interface or the constant footprint of the particles. 



Appendix A: Index contractions in complex coordinates 

In this appendix we will go into details of how to per- 
form the index contractions encountered in computing in- 
duced fields and propagator products by transforming to com- 
plex coordinates, z = x + iy and z = x — \y. We first note that 
in cartesian coordinates, covariant and contravariant elements 
of tensors are identical, i.e. the metric is the identity. There- 
fore we have not distinguished between them and have written 
all indices as subscripts in the main body of the article. How- 
ever, now that we are considering a change of coordinates, we 
will use a covariant notation instead. For example, Eq. (ITTb 
will read 



5h{r) 



a 



^., f (0)ff-"(0 ! r) 



(Al) 



where g %3 is the (inverse) metric tensor. 

From the embedding x = (l/2)(z + z),y = (l/2i)(z — z) 
the tangent vectors e z = (l/2)(l,-i) and e 2 = (l/2)(l,i) 



follow. Then the metric tensor 
found as 



e„ and its inverse are 



1 

m = o 



1 

1 o 



1 

1 



(A2) 



Since the metric is constant, the Christoffel symbols vanish, 
and therefore covariant derivatives are just partial derivatives 
(and hence also commute). 

The reason why it is useful to go to complex coordinates 
is as follows. The index contractions needed for our compu- 
tations all involve partial derivatives of the harmonic Green 
function, which is 



G{z,z') 



1 

47T 



log(z - z')(z - z ) , 



(A3) 



in complex coordinates. It is easily seen that 

d z d- z G = . 



which is equivalent to the harmonic property of the Green 
function since V 2 G = ;/' '( .', , = 4d z d z G. This is a useful 
property because, when expressed in complex coordinates, all 
index combinations involving an alternation of indices on G 
vanish from the contraction, leaving only terms that involve 
either 



d n G _ (n - 1)! 
dz n ' 



1 



4tt (z' - z) n ' 



(A5) 



or its complex conjugate. This property and the simple form 
of the metric tensor ( 1A2I ) simplify things greatly, allowing us 
to express products involving an arbitrary number of indices 
and factors. 

The computation of the induced deformation, Eq. (ITTb 
or Eq. (1AU , is now as follows: The background /i bg = 
ar cos £ip can be rewritten in complex coordinates as h bs = 
{a/2)(z e + z e ). Using Eqs. (fA^ and dA§ we find 



6h(r)=-—h^ At (0)G i ^(0,z) 



a 

1 

a 



^...r^.. l) (0)G j ,.. J( (0,z) 



2^|/i bg (0)^G(0,, 



(A6) 



Note that the index contraction, which generally has 2 e 
terms, was reduced to only two terms, owing to the prop- 
erty d z d z G = 0. Now, using Eq. dA5b the result we used 
in Eq. ( fT2l i can be found: 



Sh(r) 



CW.ail ( 1 



a 2 4tt 
(<) 2*l!(l- l)!a cosly> 

47TCT r l 



(Al) 



Note that one has to pay attention to which argument of the 
Green function is differentiated so as to avoid sign errors. 

The other point where we encounter index contractions is 
in propagator products such as 



Ql a Ql„ G abQl b ge c G bc d e c Ql d Q c<l d t d gt a q 



-id a 



(A8) 



where d e is a shorthand for d^d^ ... di c . To make the under- 
lying algebra more transparent, let us first look at this product 
of four propagators, before we generalize to an arbitrary num- 
ber. Similarly to the derivation of Eqs. dA6l l and dA7| ), and 
defining z a b = z a — Zb, this product can be written as 

m =2 t a +l M l dd i a Qt, G(Zab)d p d tc G{zbc) 

x d e /d e /G(z cd )d e /4'G(z da ) + c. c. (A9) 



(A4) Now we substitute from Eq. 



to find 
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(lA8l > 



(4 + 4 - 1)! (4 + 4 - 1)! (4 + ^ - 1)! (4 + 4 - 1)! (-2)' 

47T 4^ 4^ 4^ -A+^i. Jtb+lc Jc+td ^d+t-a ^ ^ l " 
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(4 + 4 - l)!(4 + 4 - l)!(4 + 4 - l)!(4 + 4 - 1)! cos(4^ a b e - 4v? b c d + Wa - 4^) 
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Ja +h Jb +1 c to +ld id +1 a 

'ab 'be 'cd 'da 



(A10) 



where z ab =: r ab e 1Vab and y h ac := p o6 - <^ 6c is the angle 
between r a — r b and r — r c , etc. It is obvious how this 
result generalizes to the product of an arbitrary number q of 
propagators: 

A (-2) 4 (4 + 4+i - nie^-^^-is+i 



4tt 



r 



with 4.+1 = 4- We use this formula to evaluate many inter- 
actions in Section IIIIBI most noticeably in Eq. d23l where a 
pair interaction is considered and thus the angles are zero. 

Lastly, we would like to revisit a statement we made in Sec- 
tion HHiO] namely that the free energy expansion (O for the 
surface tension Hamiltonian has identically vanishing odd cu- 
mulants. This is easily seen to hold by observing that it is not 
possible to write a nonzero expression analogous to (|A91 > for 
an odd number of propagators, due to the property (lA4t and 
the off-diagonal form of the metric. As it rests on Eq. (IA4K 
or V 2 G = 0, this property of all odd numbered cumulants 
vanishing is special to the harmonic free surface Hamiltonian. 



Appendix B: Near-contact asymptotics 

In problems where the field equation, the boundary con- 
ditions etc. are more complicated than those considered in 
this article, one generally does not expect to calculate the cu- 
mulant expansion in its entirety. However, in this specific 
case, it is possible to obtain an exact and complete expansion: 
Eq. (1231 . While finding a closed form for the whole series or 
at least the expansion coefficients is a difficult task, consider- 
able simplifications turn out to be possible. The inner 2s-fold 
sum of Eq. ( 1231 has the form of a so-called binomial cycle lf73ll 
and it is possible to perform all but one of the sums adapting 
from Ref. 17311 . Defining u = R 2 jr 2 for the case of identical 
radii, Eq. d23l reduces to 



'2n\ u 2s+n+2 7"(*-n,l; 



s,n— 



72 



' hs+i(u) 



(s + 1)/ 2s +i(u)/ 2s (m) 



(Bl) 

where J 7 (a, b; c; x) is the hypergeometric function 117411 and 
the polynomials f q (u) have the recursion relation 



fq(u) = fq-l(u) - U/ 9 _ 2 (U) , 



(B2) 



with /o = /_i = 1, related to continued fractions. Also note 
that the summation index s merely labels cumulants and is 
related to the cumulant order q as q = 2(s + 1). This already 



amounts to a big practical improvement over Eq. d23l , because 
Eq. (IB lb can be much more efficiently expanded in powers of 
inverse distance. 

From the recursion relation of f q , one can also derive the 
expression 



/«(«) = E 



i=0 



q+l-i 



(-uY , (B3) 



where [_-J denotes the integer part, or better still, recognize (in 
our case, MATHEMATICA did this for us) its closed form as 



/,(«) 



X 



9+2 _ l 



(x 2 -i)(x 2 + i) 9 ' 



(B4) 



where x(u) = (1 — y/1 — Am — 2u)/2u. Thanks to the closed- 
form expression for f q , one can show that in the interval 
< u < 1/4 (i.e. from asympotically separated to osculating 
disks, recalling u = R 2 /r 2 ), the ratio /2s(w)//2s+i(w) < 2, 
which is the condition for the hypergeometric function in 
Eq. ( IBU to converge as n — > oo: 



(B5) 



One can then use Eq. ( IB5I ) to isolate the leading divergence in 
each cumulant q = 2(s + 1) in Eq. ( IB II ) as the n — > oo tail of 
the inner sum and obtain the resulting divergence in j3F as 



lim J r (— 71, 1; — 2n; x) = I 1 — 



/3F a 



OO 

■E 



.2s+2 



n=0 \n) 



(B6) 



The symbol "~" reminds us that this only captures the leading 
order asymptotic divergence. 



After noting that J^l 



(2n\ 



= (1 — 4u)~2, we will 
define the (scaled) surface-to-surface separation d between the 
disks as u = l/(d + 2) 2 and re-express Eq. JB6I > for d< 1 
using Eq. ( IB4I ) in order to obtain the leading divergence in 



OO - 

72 



C(2) _ 7T 2 

4Vd 24Vrf 



(B7) 



in agreement with Eq. 

The crudeness of isolating the divergent part only allows 
one to obtain the leading order proximity asymptotics. In a 
recent preprint ll70ll one of us used conformal field theory to 
exactly compute the partition function of a bosonic field on a 
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plane with two holes. This is isomorphic to the problem of 
two disks on a film studied in this article, and the interaction 
can in fact be written in closed form (it involves a Dedekind 
jy-function). Its expansion at large distances reproduces our 
EFT series, while an expansion at contact leads to Eq. ( |25| >. 

We again point out that computing the close distance diver- 
gence is not the most obvious application for our specific EFT 
implementation (which is built on a large distance expansion), 
but it nevertheless illustrates the power of the formalism quite 
vividly. 



Appendix C: Cancellation of divergences by counter-terms 

In this appendix we will illustrate how self-energy diver- 
gences can be removed by counter-terms. While this is, to 
some extent, a standard textbook affair, we include it in order 
to illustrate the mechanism in the presently relevant context 
for readers not necessarily familiar with (quantum) field the- 
ory. 

First of all, let us clarify which divergences are of inter- 
est here. Purely numerical divergences, such as those of the 
first cumulant, get absorbed in the definition of the free sur- 
face energy (i.e. no particles); they do not refer to the relative 
positions of the particles and thus do not affect interaction en- 
ergies. Similar trivial divergences occur in every cumulant 
order and do not require any further discussion. The more 
interesting divergences are those which multiply interactions 
between particles and thus cannot be lumped together with the 
free surface energy. 

The appearance of these divergences is due to the very con- 
struction of an effective field theory. The self-interactions 
necessarily involve the short distance (high energy) physics, 
through G(0), where the theory is inadequate; especially in 
field theories of condensed matter, we know that the contin- 
uum description breaks down below a lattice-spacing or an 
equivalent short length scale. But the effect of the short dis- 
tance physics on large length scale (low energy) observables, 
such as the polarizabilities is to renormalize (or "dress") 
their values. The values we obtained for these polarizabilities, 
through a procedure that has nothing to do with the short dis- 
tance physics of the field, are therefore renormalized values. 
On the other hand, if a field theory is to make finite predic- 
tions on (low energy) physical observables, the renormalized 
couplings must be accompanied by counter-terms, or equiva- 
lently the coupling constants must be restored to their "bare" 
values. However, since the counter-terms exist for the sole 
purpose of canceling these ultraviolet divergences in the val- 
ues for physical observables, it is legitimate not to explicitly 
write the counter-terms, but discard such divergences on these 
grounds when encountered, as we did in this article. We will 
nevertheless provide an example below. 

For the sake of simplicity, let us assume two particles pos- 
sessing only a dipole polarizability described by the Hamilto- 
nian 



ah 



+ l -C 2 hf{r 2 ) 



(CI) 



It is not hard to see that the divergences in the first and sec- 
ond cumulants, (PA'H) C and -(1/2) ((PAH) 2 ) c , are of the 
trivial sort, i.e. those that do not involve the distance between 
particles. However, in the third cumulant we encounter a di- 
vergent pair free energy 



1 

3! 



G 2 G 2 
(2a)3 



Q/"~vll /nfl2 s~i'21 



(C2) 



as well as a similar term with the labels 1 and 2 interchanged. 
Using the fact that ^5^(0) = (1/2)<%V 2 G(0), we may 
rewrite (IC2t as 



^jk^kj 

(7 J 



(C3) 



where it is apparent that this energy scales as (5(0) /r 4 . The 
way the factors were arranged in dC3l ) makes it easy to iden- 
tify what counter-term should be added to the Hamiltonian to 
negate this unphysical divergence: One might check that the 
physical pair energy between these two particles at the second 
cumulant has exactly the form of ( IC3b . except that the factor 
in square brackets is C\ja. This means that if one adds to 
ATi the counter-term 



1C? 

2 2(7 



V 2 G(0)/i?(n 



(C4) 



it produces a divergent term in the second cumulant that ex- 
actly cancels the aforementioned divergence, ( IC2b or ( IC3b . in 
the third cumulant. In fact, this counter-term removes all di- 
vergences caused by a single occurrence of G\j in the entire 
cumulant expansion , but it is beyond the scope of this article 
to prove this. Also note that adding this counter-term to the 
Hamiltonian is equivalent to redefining the dipole polarizabil- 
ity C\ as 



Ci -> C x 1 



^V 2 G(0) 

Z(T 



(C5) 



This was a simple example that exhibits an induced dipole 
charge that produces a divergence due to interacting with it- 
self once. Considering other multipole orders will introduce 
new counter-terms involving different divergences, such as 
V 4 G(0). Removing divergences due to one charge interact- 
ing with itself more than once is achieved by counter-terms of 
higher orders in the divergences, such as (V 2 G(0)) 2 . Regard- 
ing these counter-terms as restoring the polarizabilities to their 
(divergent) bare values C^L from the renormalized one 
may write 



_ r {n) 

U bare — ° 



/ j n 
n=l 



where the coefficients b n are pure numbers appropriately 
chosen to cancel the offending diagrams. In this appendix, 
for instance, we have explicitly found out that = 1/2. 



V^G(O) 



(C6) 
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Appendix D: Induced monopoles 

In our treatment we avoided induced monopole terms in the 
effective theory, i.e. terms of the form 



(Dl) 



since the field of these charges — each proportional to log|i — 
r a \ — violates square integrability of its gradient in M 2 . Such 
terms would describe particles with frozen vertical fluctua- 
tions. We will discuss a possible workaround for this issue. 
Consider the regularized free surface Hamiltonian 

H[h] = \° j d2r { h l + X~ 2 h 2 ) (D2) 

which approaches the tension Hamiltonian (ffj in the limit of 
large A. For an interface between two fluids subject to gravity 
(horizontal in its unperturbed state), A is the capillary length 
and is given by A = dajgp, where g and p are the gravi- 
tational acceleration and mass density difference between the 
fluids, respectively. The addition of the "mass" term damps 
correlations over distances larger than A, hence regularizing 
the infrared divergence of monopole fields. The Green func- 
tion for this choice of surface energy is 



G(r,r') = -Liifo(|r-r'|/A) 



(D3) 



instead of ( — l/27r) log \ r — r'\, where Kq(x) is a modified 
Bessel function of the second kind. After this regularization, 
we may safely consider particles that are completely pinned. 

The effective theory of such particles involves the same in- 
duced charges as before, Eq. (TT3T >. with the addition of the 
monopole (£ = 0) terms dD lb . The reader could object that, 
with the new choice of kernel, there may be new terms such as 
/i 2 or hijkkhij, etc., but these terms can be eliminated using 
the equation of motion, (—V 2 + X~ 2 )h = (for a proof, see 
e.g. Ref. 1I57I0 . The matching of the polarizability coefficients 
is done similarly by comparing the response of the induced 
charges to backgrounds of the form /i bg = ar cos £<p in the 
full and effective theories. One finds 



5h tull {r) 



ft 



Ki(R/X) 



Ke(r/X) cos £ip 



and 



Sfi 



off I 



C £ £\ 



a 2tt) 



■Ki(r/X) cosif , 



respectively, yielding 

C (l) 



£\K e (R/X) ' 



(D4) 



(D5) 



(D6) 



Expanding the Bessel functions, one observes that for £ > 
these polarizabilities converge to those in Eq. ( fT4t as A —> 00, 
whereas the massless limit of the monopole polarizability is 



c (o) 



2na 



e 1 R 



as A 



(D7) 



where 7 is the Euler-Mascheroni constant. 

We observe that the monopole polarizability vanishes like 
1/ log(i?/A) in the massless limit of infinite capillary length 
that we eventually want to take. This means diagrams in- 
volving monopoles could vanish as well, unless every factor 
of 1/ log(i?/A) due to a monopole polarizability is canceled 
by a similar factor in the numerator. Factors of this form 
indeed exist: they come from monopole-monopole links in 
the propagator product, since in the massless limit G(r) = 
— (l/27r) log(e 7 r/2A). Due to the closed topology of the di- 
agrams, there are enough monopole-monopole links to bal- 
ance the vanishing monopole polarizabilities only when there 
are no higher order multipoles in the diagram; replacing one 
monopole polarizability in the diagram costs two monopole- 
monopole links. In other words, the only monopole inter- 
actions that do not vanish in the massless limit are those 
with other induced monopoles and nothing else. This elu- 
cidates and generalizes the findings of Lehle and Oettel that 
monopole-dipole and monopole-quadrupole interactions in- 
deed vanish M45I1 . 

We can now write the pair interaction between two pinned 
particles, on a capillary surface for which the capillary length 
tends to infinity. The interaction will consist of Eq. d23l l due 
to induced multipoles of order £ > and, based on the dis- 
cussion of the previous paragraph, a part f3F mon that is due 
solely to monopole polarizabilities. Observe that in the latter, 
propagators do not carry any derivatives and therefore all the 
cumulants are of the same order in the inter-particle separa- 
tion. Hence, all cumulants must be summed for the monopole 
interactions. To evaluate this, one can refer to Eq. ( 1221 . keep- 
ing in mind that there is no sum over multipole orders £i now; 
only the monopoles are taken. Assuming identical particle 
radii R to declutter expressions, one finds 



/3F D 



lAl f C^G{r) 
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s \Ko(R/X) 
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Owing to the damping of the regularized theory, this series 
converges for all r (<C A — > 00) to 
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lop 



1 - 



](r/X) 



-— l0E 

2 e 
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K Q (R/X)-K (r/X)} 

K Q (R/X) + K Q (r/X)} 



const. 



(D9) 



When the massless limit A — > 00 is taken, the first term on the 
right hand side of the last equality gives 



- log log , 



(D10) 



in agreement with Ref. [14511 . The second term is propor- 
tional to loglog(A 2 /i?r), which is associated with a force 
~ 1/r \og(Rr / X 2 ) — s- as A —> 00. We note that extension to 
particles of unequal radii results in the change R y/RiR 2 
in Eq. dDlOl 
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